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ABSTRACT 

We report the discovery of a large timescale candidate microlensing event of a bulge stellar source 
based on near-infrared observations with the VISTA Variables in the Via Lactea Survey (VVV). The 
new microlensing event is projected only 3.5 aremin away from the center of the globular cluster 
NGC 6553. The source appears to be a bulge giant star with magnitude = 13.52, based on the 
position in the color-magnitude diagram. The foreground lens may be located in the globular cluster, 
which has well-known parameters such as distance and proper motions. If the lens is a cluster member, 
we can directly estimate its mass simply following iPaczvnski et ahl (1199611 which is a modified version 
of the more general case due to Refsdal. In that case, the lens would be a massive stellar remnant, 
with M = 1.5 —3.5M0. If the blending fraction of the microlensing event appears to be small, and this 
lens would represent a good isolated black hole (BH) candidate, that would be the oldest BH known. 
Alternative explanations (with a larger blending fraction) also point to a massive stellar remnant if 
the lens is located in the Galactic disk and does not belong to the globular cluster. 

Subject headings: globular clusters: individual (NGC 6553) - gravitational lensing: micro - infrared: 
stars - surveys 


1. INTRODUCTION 


Microlensing is the best (and so far the only) tool 
to detect isolat ed stellar mass black holes (BHs; see 
iMao et ^ l20I2tl . The microlensing projects have dis¬ 
covered a few can didate stellar-mass BH lenses toward 
the Galactic bulge (iBenne tt et al.[ )2002) :lMao et al.ll2012l: 
lAgol et aI1l2002l : TPoindexter et ahKOOSll . I Gould fc SalimI 
( I999f) sp^ificallv applied astrometric microlensing in or¬ 
der to measure the masses of BHs. Although the mea¬ 
surements remain challenging, this method is importan t 
fo r future measu rements from space (|Gould fc Yeell201^ . 

iRefsdall ()1964[ 1 showed that the mass M of a gravita¬ 
tional lens is related to observables {prei,i^rei,tE) by: 
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where Ie is the Einstein crossing time derived from 
the event light curve and Ttrei and prei stre the lens- 
source reUtiw parall ax and proper motion respectively. 
While IRefsdall (jl964[ l clearly ha d in mind that t h e lens 
and source could both be seen, IPaczvnski et al.1 (119941 1 
pointed out that for globular-cluster lenses detected via 
their magnification of Galactic-bar sources, iTrei and prei 
can potentially be estimated without directly observing 
the lens. This opens the possibility that the method 
could be applied to BH lenses. The main advantage of 
such a search is that the microlensing geometry is more 
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certain, and there is a direct relationship between the 
lens mass and the time scale of the event. This tech¬ 
nique was successful in the case of the first confirmed 
globular cluster microle nsing event in M22 (NGC 6656 ) 
due to a low-mass star (|Pietrukowicz et aI1l2005l [201211 . 
Later, another microlensing event in close pro jection to 
the g lobular cluster NGC 6553 was reported (iVee et al.1 
I2013fl . which also turned out to be a low-mass lens, 
possibly with a plane t H owever, the suggestion by 
IPaczvnski et al.l (11994 1199611 seems to be a promising 
way to search and measure heavy remnants (neutron 
stars (NSs), or BHs) in spite of the technical challenges. 

Here we report the discovery of a long-timescale can¬ 
didate microlensing event of a bulge stellar source pro¬ 
jected only 3.5 aremin away from the center of the glob¬ 
ular cluster NGC 6553. If lens membership to the cluster 
is confirmed, the lens could be a massive stellar remnant, 
possibly a BH with ~ M = 2Mq. Section 2 presents the 
VVV Survey observations of the microlensing event. Sec¬ 
tion 3 discusses the globular cluster NGC 6553 where the 
lens may reside. Section 4 shows the measured parame¬ 
ters of the lens, and Section 5 summarises the conclusions 
and future tests. 

2. VVV SURVEY OBSERVATIONS OF THE 
MICROLENSING EVENT 

VISTA Variables in the Vfa Lactea (VVV) is a public 
ESO near-infrared (IR) var iability survey aime d at scan¬ 
ning the inner Milky Way (iMinniti et al.l[2M^ . The ob¬ 
servations are acquired with the VISTA 4 m telescope 
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TABLE 1 

VVV Survey near-IR photometry of the candidate 
MICROLENSING EVENT. 


Epoch (MJD) 

Ks 


56363.376890 

13.509 

0.013 

55298.333635 

13.509 

0.015 

55491.015648 

13.528 

0.013 

55778.139393 

13.495 

0.013 

55804.103650 

13.543 

0.012 


Note. — Table 1 is published in its entirety in the electronic 
edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

at ESO Paranal Observatory ([Sutherland 1 1201^ . The 
VVV survey covers an area of 562 sqdeg in the inner disk 
and bulge of the Milky Way. The VVV database now 
contains multicolor [ZYJHKs) photometry, and multi¬ 
ple epochs in th e JCs-band, moni toring about a billion 
sources in total (|Saito et al.ll2012D . The i^s-band obser¬ 
vations continue, and the variability light curves so far 
span from 2010 to 2015. This large database enables a 
number of studies of different variable ob jects, including 
micro lensing events ([Catelan et al.l I2013t iHemoel et al.! 

[Ml . 

We have started to search for variable stars using 
the VVV Survey data in the fields of 36 globular 
clusters toward the Galactic bulge. The data reduc- 
tions and photometry have been described in detail by 
lAlonso-Garcia et al.l (j2015li . where our first results for 
the globular clusters 2MASS-GG02 and Terzan 10 are 
also presented. Briefly, the individual VVV Survey im¬ 
ages are reduced, astrometrized and stacked by the Gam- 
bridge Astronomy Survey Unit (GASU) u sing the VISTA 
Data Flow System (yDFS) pipeline ([Emerson et al.l 
12004 llrwin et M l20f)4 iHamblv et al.l 12001 . and the 
photometry is calibrated onto the VISTA filter system. 
Then we carry out point-spread function (PSF) photom¬ 
etry on the individual processed stacked images using 
an updated version of Do PHOT (jSchechter et al.l 119931 : 
lAlonAGarcia et al.ll2012li . The photometry for each ob¬ 
ject in the different images available was cross- correlated 
using the STILTS package [Tavlor et al.l (|2n06l i. and the 
ligh t curves generated were the n analyzed for variability 
(see lAlonso-Garcia et 3111201511 . 

During this near-lR search for variable stars 
in the field of the glob ular cluster NGC6553 
([Gontreras Ramos et al.ll2015l in preparation), we dis¬ 
covered a candidate microlensing event located at 
R.A.(2000)=18:09:13.86, decl.(2000)= -25:57:52.7, that 
is only 3.5 arcmin away from the cluster center. This 
event peaked in the 2012 bulge season (on the 2012, July 
8), and it has not been reported by any of the ongoing 
microlensing surveys. The brightness of the object in¬ 
creased by about 0.7 mag in the iLg-band over 50 days, 
reaching a magnitude of Ks = 12.8, and then fading to a 
constant level of ATg = 13.5. Table 1 shows the ATg-band 
observations for this object. The ATg-band epochs corre¬ 
spond to different nights, and the seeing for the majority 
of them was very good, between 0.6 and 1.2 arcsec. 

Fig. 1 shows the 6.8 x 3.5 arcmin finding chart of the 
microlensing event with respect to NGG6553. This lies 
well within the cluster tidal radius of A = 8.16 arcmin. 
The zoomed 30 x 30 arcsec^ region to the right shows that 
the source is located just in between two brighter stars 


(both saturated with ATg < 12) that are separated by 3 
arcsec. This makes the photometry difficult in the worst- 
seeing images, explaining a few of the deviant points with 
large error bars in the light curve (Fig. 2). 

The source appears unblended in the best seeing im¬ 
ages (we have also checked in the Hubble Space Tele¬ 
scope (HST) archive and unfortunately there is no im¬ 
age of this field), and Fig. 2 shows the light curve for 
this eve nt, along with the simple un blended microlens¬ 
ing fit (jPaczvnski et al.l [l994 TTo^ . This would cor¬ 
respond to an event without significant contamination 
from an unresolved source. However, we have fitted a 
microlensing light curve, obtaining two equally good fits 
with = 200: an unblended and a blended fit, that are 
discussed in turn in the following sections. Both fits are 
very good, with the following microlensing parameters: 
Gase A (without blending): baseline magnitude = 
13.515 ± 0.002, impact parameter ug = 0.62 ± 0.004, 
time of closest approach tg = 56117.5 ± 0.43, Einstein 
timescale ts = 51.3 ± 0.8 days, blending / = 1.00, and 
= 201. Gase B (with blending): baseline magnitude 
Ks,o = 13.515±0.002, impact parameter ug = 0.46±0.08, 
time of closest approach tg = 56117.4 ± 0.42, Einstein 
timescale ts = 62.5 ± 9 days, blending / = 0.61 ± 0.17, 
and = 199. 

The values obtained for show that there is no im¬ 
provement in the curve fitting when a blending factor is 
included. Also, the VVV Survey is designed to monitor 
variable stars, and the light curve is not frequently and 
evenly sampled, preventing us from using more complex 
models (with larger numbers of free parameters like those 
including the parallax effect) with better confidence than 
a simple microlensing light curve. 

Can this be a variable star instead of a microlensing 
event? The main kinds of variable star contaminants are 
foreground dwarf novae and distant background SN, both 
of which are unlikely. The shape of the light curve is 
very symmetric, unlike the nova, dwarf novae or super¬ 
nova light curves that typically show a fast rise and a 
slow decline. There are some exceptions, as some re¬ 
current novae show symmetric outbursts (see for exam¬ 
ple, the light curve of GK PeiQ). However, we discard 
a dwarf nova outburst because of the light curve shape 
and also because our object is too red (J — ATg = 1.0) 
to be a dwarf novae. We also searched for variability in 
the region outside the peak but no signal of eclipses or 
ellipsoidal variations from a binary system were found as 
expected from a cataclysmic variable. Also, even though 
the VVV Survey data is deep enough to see background 
galaxies throughout the bulge, in this case no extended 
galaxy-like object that could be a supernova host galaxy 
is detected. 

3. THE GLOBULAR CLUSTER NGC 6553 

Because of its proximity to a globular cluster, this 
microlensing event is particularly interesting. The 
bulge globular cluster NGG6553 has been well studied, 
and its physical parameters are relatively well known. 
NGC 6553 is an old a n d metal-rich bulge g l obular cluster 


(lOrtolani et IMinniti et alJll99r^ 1 

Barbuv et alJ 

119981: IZoccali et al.ll2001l lAlves Brito et al. 

1200611. This 


cluster is moderately concentrated with a core radius 


^ http://www.aavso.org/sites/default/files/images/LTGKper.gif 
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Fig. 1.— Left: 6.8 x 3.5 arcmin finding chart showing the position of the microlensing event with respect to NGC6553. The projected 
distance to the cluster center is 3.5 arcmin, well within the cluster tidal radius oi R = 8.16 arcmin. Right: zoomed 30 x 30 arcsec^ region 
centered on the source star, showing that it is located just in between two brighter stars (with Ks < 12) that are separated by 3 arcsec. 
The seeing of the images is 0.8'^ 



Time (MJD) 

Fig. 2.— i^s-b and light curve for this event, along with a simple 
microlensing fit (|Paczvnski et a0ll994l') . The microlensing param¬ 
eters of this fit are: baseline magnitude Ks,o = 13.515 zb 0.002, im¬ 
pact uq = 0.62di0.004, time of closest approach to = 56117.5di0.43, 
Einstein timescale = 51.3 ±0.8, blending fraction / = 1.00, and 
= 200. The points with the larger errorbars correspond to the 
worst seeing images. 

Tc = 0.55 arcmin, and a tidal radius rt = 8.16 arcmin. 
It is very reddene d, with E(B — V) = 0.73, Ay = 2.26, 
and Ak„ = 0.23 l|Barbuv et aLlHO^ . Particularly rel¬ 
evant for this work is the distance of the cluster. We 
adopt a distance of 6.0 kpc for this cluster following 
I Alves Br ito et a 1 (1200611 . as listed in the 2010 version 
of Harris et al.l ( 199611 catalog, noting that the distance 
uncertainty adds to the total error budget. We observe 
the cluster horizontal branch (HB) at = 12.5 (Fig. 3), 
w hich is consistent w ith this distance measurement. 

IZoccali et al.l (120011) measured a relative proper motion 
of NGC 6553 with respect to the bulge of fi£ = 5.89 and 
/ih = 0.42 mas/yr. This gives a relative mean proper 
motion difference between bulge and cluster stars in the 
sky of 5.9 mas/yr. That would be the expected mean 
relative proper motion of the bulge source and the lens 
if the latter is a globular cluster member. Of course, the 
large bulge velocity dispersion can result in a larger or 
smaller actual relative proper motion. 

The source star in particular (with = 13.5, J — 
Ks = 1.0 at the baseline) is redder than the RGB of 
the globular cluster, consistent with a bulge giant. The 
source is also fainter than the bulk of the bulge red clump 
stars. If the source is a red clump giant star, it must be 



Fig. 3.— Ks. vs J — Ks near-infrared color-magnitude diagram 
(CMD) of about 20000 stars in 4 x 4 arcmin^ field centered on the 
globular cluster NGC 6553. The red circle marks the source star 
of the microlensing event. This CMD shows the foreground disk 
main sequence, the populated globular cluster main-sequence (MS) 
turn-off, the cluster red giant branch (RGB), its red horizontal 
branch (HB), and its red giant bump, as well as a redder and wider 
background bulge RGB, including the bulge red clump. We caution 
that photometric non-linearity and saturation starts at Ks < 12.5. 

located on the far side of the bulge, at Dg ~ 9 kpc. 

The field is very crowded and contains numerous bulge 
giants, and therefore membership of the lens to the clus¬ 
ter cannot be unambiguously secured. Luckily the VVV 
Survey images cover a large area (1.5 sqdeg), and we 
can measure the background contamination. At the dis¬ 
tance of 3.5 arcmin from the center of the cluster, we 
estimate that roughly 20% of the stars per unit area are 
cluster members. However, that may or may not apply 
to the distribution of heavy re mnants. In contr ast, the 
microlensing event discussed bv lYee et al.l (j2013D was lo¬ 
cated 6.7 arcmin away from this cluster, where the clus¬ 
ter stellar density is 6%, based on which they estimated 
that there is only 18% probability that their lens is a 
cluster member. Following their procedure we estimate 
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that there is roughly 50% probability that the lens, lo¬ 
cated 3.5 arcmin from NGC 6553, actually belongs to the 
cluster. 

4. THE CANDIDATE BH LENS 

The dark microlensing events toward the bulge (i.e. 
with no light from the lens) that have long-duration 
(> 100 days) may be due to massive stellar remnants. 
However, the masses based on timescales are statistical 
mass estimates. Long-duration events can also be due to 
slow-moving, low-mass objects. This makes the present 
event particularly interesting, since an Einstein diame¬ 
ter crossing time of t = 103 days (jBennett et al.l (j2002lB . 
may be due to a relatively massive lens, such as a heavy 
non-luminous remnant (NS or BH). 

Because of the position of the source star in the CMD 
(Fig. 3) we can assume that it is located in the Galactic 
bulge. Under this assumption, we consider two scenarios: 
In the first scenario the lens is a member of the globular 
cluster N GC 6553. In this case we can estimate the lens 
mass fc.f. iPaczvnski et 31111994 Il996ll . because we know 
the distance and the proper motion of the cluster. 

We have fitted different possibilities considering differ¬ 
ent distances to the source star. The mass determination 
is uncertain due to the unconstrained Urp.i ■. due also to th e 
considerable bulge line of sight depth (iNataf et al.ll2013ll . 
The most likely source distance is where the density of 
bulge stars peaks along the line of sight, i.e. Dg = 8 kpc. 
In this case, assuming Dl = 6 for the lens, we obtain a 
lens mass Ml = 3.5 ± 0.1 Mq (using relative transverse 
velocity 220 km/s), or Ml = 2.0 ± O.IM 0 (using rel¬ 
ative transverse velocity 168 km/s). If the source star 
is located in the far side of the bulge, at Dg = 9 kpc, 
we obtain a smaller lens mass, Ml = 2.7 ± 0.1 Mq and 
Ml = 1.5 ± 0.1 Mq for both cases, but still within the 
realm of the heavy remnant hypothesis. As discussed 
in Section 2, there are two equally good fits, one with 
a dark lens (case A with blending parameter / = 1.0), 
and another one with a significant contribution from an 
unresolved source (case B with / = 0.62 ± 0.20). In the 
case of the first fit, if the lens is a member of NGG 6553, 
we estimate M = 3.1 ± 0.4Mq and M = 1.8 ± 0.2Mq for 
both cases, where the uncertainty is driven by the source 
distance and transverse motions. This fit then rules out 
other stellar remnants such as white dwarfs, and is more 
consistent with a NS or BH stellar remnant. 

Even though we claim priority for recognising this as 
a BH candidate, this event was independently found by 
OGLE0as event OGLE-2012-BLG-0548. We would like 
to note that the reported OGLE microlensing parameters 
agree quite well with the ones derived solely from VVV 
p hotometry. 

IBennett et al.l (1200211 estimate that the mean mass for 
six microlensing parallax events that they study is 2.7 
Mq, arguing that they are BH candidates because they 
surpass the Chandrasekhar mass of 1.4 Mq, being even 
larger than the allowed maximum NS mass of 2.0 Mq. 
Interestingly, the mass we measure for the lens studied 
here is similar to the mean mass e stimated for the six lo ng 
timescale microlenses studied bv IBennett et al.l (| 2002 ll . 

This is an interesting case to study, because if the lens 

^ reported online at: 

http: / / ogle.astrouw.edu.pl / ogle4/ews/2012/ews.html 


is SO massive and old, it should be located close to the 
center of the cluster. It is unexpected to find a heavy 
stellar remnant at a distance of about 6 core radii from 
the cluster center so long after its formation. For exam¬ 
ple, the populations of millisecond pulsars and LMXBs 
in globular clusters tend to be among the most radially 
concentrated, and it is not clear why this BH did not 
sink into the globular cluster ce nter through dynamical 
friction ('e.g.. iHeinke et al1l2005ll . 

As another point, we note that XMM source #27 
in NGC 6553 from IGuillot et al.l (120 III! located at 
R.A.(2000)=I8:09:2I.67, decl.(2000)= -25:57:31.6 is the 
closest X-ray source to the position of the microlensing 
event located about 2 arcmin away. Based on the existing 
data, however, it appears unlikely that XMM source #27 
corresponds to radiation from this candidate BH, as ex¬ 
pected if it is accreting gas from the i nterstellar medium 
or the wind of a companion star fe.g.. lAgol et al.l[2002ll . 

In the second scenario that we consider, the lens is not 
a member of the globular cluster NGC 6553. In this case 
the distance to the lens is unknown and its mass is not 
well constrain at all. However we can make some as¬ 
sumptions to visualize possible solutions. If we consider 
the distance of the bulge for the source, a long timescale 
can be obtained for a low-mass star as lens only if this 
lens is located at a few kpc from the Sun. For a source 
at a distance of Ds = 8 kpc, and assuming that the lens 
is located at Dl = 3, 4, and 5 kpc with a typical disk 
velocity of U = 220 km/seg we estimate: M 3 = 2.8Mq, 
M 4 = 2.6M©, Ms = 2.8Mq, respectively. Instead, for 
a source at a distance of Ds = 9 kpc, and a lens at 
Dl = 3, 4 kpc, and 5 kpc we estimate: M 3 = 2.6Mq, 
M 4 = 2.3Mq, Ms = 2.3Mq, respectively. 

This case would favor more the second fit, with signif¬ 
icant blending from an unresolved source. However, this 
is interesting because the inferred masses are large for 
a typical disk main sequence star. If the lens is a main 
sequence star, it should be bright enough for detection 
at these intermediate distances, and also bluer than the 
observed color J — Kg = 1.0, making this scenario less 
likely. A lens closer to the Sun than 3 kpc is unlikely be¬ 
cause of the same argument: it would be bright enough 
and should have been observed. 

5. CONCLUSIONS 

We have re-discovered a microlensing event in the field 
of the globular cluster NGC 6553. Because the distance 
of this cluster is well known {D = 6.0 kpc), and the 
CMD suggests that the source is likely a bulge red gi¬ 
ant at Ds = 8 — 9 kpc, we can solve for the mass of 
the leasing object if we assume this is a cluste r member. 
Simply using eq. 15 of IPaczvnski et all (1199611 results in 
a lens mass of M = 1.5 — 3.5Mq consistent with a heavy 
remnant. 

There is another equally good fit to the microlensing 
light curve, with significant blending from an unresolved 
source. However, we note that combining the optical 
and IR data from OGLE and VVV, the microlensing 
color can be more precisely measured, and it appears 
to be essentially similar to the baseline color, meaning 
that the blending is really small, or alternatively that 
the lens has exactly the same color as the source. Proper 
motion measurements are warranted in this case, and 
combining VVV and OGLE data would greatly increase 
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the precision of the parallax measurement and will be 
of great importance to better constrain the mass of this 
object. The alternative scenario of a more nearby disk 
main sequence star leasing a distant bulge giant cannot 
be ruled out, but also yields a rather massive lens (2.3-2.8 
Mq). Even though we cannot rule out this possibility, 
it appears to be less likely because in this case the lens 
would be detectable as a bright blue star. 

We can make specific predictions if the leasing object 
is a cluster member: the lens should eventually move 
away from the source because the relative proper motion 
of this globular cluster is relatively large, and then be 
detected if it is a stellar object, or remain invisible if it is 
an isolated BH. This method is the best way to confirm 
the m icrolensing nature of the event iPietrukowicz et al.l 
(|2012ll . In the present case, the source-le ns separation 
can r each the 60 mas range within 10 year (IZoccali et al.l 
1200111 , and the leasing star should be detected with high- 
resolution images (with HST or AO cameras) if it is not 
a BH in the cluster. In this sence, immediate high reso¬ 
lution A'-band measurements using AO are desirable in 
order to provide a reference image for subsequent detec¬ 
tion of the lens. 

If membership to the cluster is confirmed with follow¬ 
up observations, this would not only be the only known 
stellar mass BH in a globular cluster, an important ob¬ 
ject to validate stellar evolution theory, but also the old¬ 
est BH known. This discovery opens up some interesting 
questions, like what is the mass distribution of globu¬ 
lar cluster BHs produced by the death of ancient mas¬ 
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